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1 Introduction and Results 

One of the interesting experiments in particle physics is the heavy ion colli- 
sion done at RHIC (Relativistic Heavy Ion Collider) and LHC (Large Hadron 
Collider). In these experimental set-ups two pancakes of heavy nuclei such 
as Gold(Au) or Lead(Pb) are collided at a relativistic speed. Very soon af- 
ter the collision a new phase of matter called QGP (Quark-Gluon Plasma) 
is produced. Understanding the properties of the QGP has attracted a lot 
of attention. It's been realized through hydrodynamic simulations and ex- 
perimental observations that the plasma is strongly coupled with very low 
viscosity over entropy density, (^) [1,2]. 

A very challenging observation is the very rapid thermalization of plasma. 
At almost Ifm/c after the collision the plasma reaches local equilibrium (ther- 
malizes) and its dynamics is approximately described by ideal fluid hydrody- 
namics [3,4]. Since the system is strongly coupled and perturbative calcula- 



tions can not be applied, it is difficult to study an out-of-equilibrium process 
such as thermalization. Models of QCD such as lattice QCD have not been 
able to achieve much to explain this phenomenon. 

Recently AdS/CFT correspondence [5], or more broadly, the gauge/gravity 
duality has gained some success in making toy models to describe this pro- 
cess. Generally speaking, the AdS/CFT duality has provided a powerful 
framework to study strongly coupled gauge theories. According to the origi- 
nal statement of the correspondence, type JIB string theory on AdS^ x S^ is 
dual to D = A, Af = A, SU{N) super Yang-Milles (SYM) theory. This idea 
was then generalized to a broader class of gauge/gravity dualities in various 
dimensions [6]. It was shown that in the large A^ limit, the strongly coupled 
SU{N) gauge theory living on the worldvolume of A^ coincident Dp-branes 
is dual to the supergravity on the near horizon geometry of the Dp-branes. 
Although for p = 3 the gauge theory enjoys conformal symmetry, the gauge 
coupling is dimensionful for the other values of p. 

The absence of conformal invariance in the gauge theory leads to the 
radial variation of both the string coupling (dilaton) and the spacetime cur- 
vature in the dual description. This is due to the fact that duality relates the 
radial coordinate transverse to the Dp-brane to the energy scale in the gauge 
theory. The supergravity background is only reliable for weak string coupling 
and small curvature. This is provided in a regime of energy where [6] 

1 <^e// < A^^ . (1.1) 

geff is the dimensionless effective coupling, g^ff = gyj^jNu^''^, where u is 
the radial coordinate. In this intermediate regime of energy the dual gauge 
theory is always strongly coupled. 

AdS/CFT techniques have been applied to study rapid thermalization ob- 
served at RHIC [7-11]. The thermalization in field theory which is an out-of- 
equilibrium process happens after an injection of energy into the system. In 
the gauge/gravity context this injection can be done in two different ways: by 
directly adding a time-dependent source to the boundary field theory which 
corresponds to turning on a non-normalizable mode [7,8] or by introducing a 
time-dependent coupling in the field theory which is represented by nontriv- 
ial time-dependent classical solutions of some kinds of probe branes [10, 11]. 
In the ffist method the source is assumed to be nonzero only for a limited 
time interval. In the bulk gravity dual this corresponds to the collapse of 
matter and formation of the horizon; in other words hlack hole formation. 



In [7] the authors have been able to quantitatively give an approximate value 
for the isotropization (thermalization) time of the plasma which is close to 
the experimental result. In the other method the horizon formation happens 
on the probe brane and then the energy can flow from the probe brane to the 
bulk gravitational degrees of freedom. This corresponds to the dissipation of 
energy into the field theory. 

Note that in the first method the field theory is pure YM and conformal. 
While in the other one the field theory contains fundamental matters and 
mesons in addition to the pure YM. Studying the thermalization in the me- 
son sector is the subject of an interesting work [11] where the authors have 
modeled the thermalization by baryon injection into the system. This idea 
comes from the observed sudden change in the baryon chemical potential at 
QGP production. 

The addition of fundamental matter to the SYM theory which lives on 
the boundary of AdS space is done by adding the fiavour branes to the dual 
gravity system. This has been proposed for the first time in [12] where Nf 
D7-branes are added to the AdSs x S^ background in the probe limit where 

is small. This will reduce the number of supersymmetries of the original 
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theory from 16 to 8 and will add a matter hypermultiplet in the fundamental 
representation to the SYM theory. The hypermultiplet describes the dy- 
namical quarks living in four dimensions and the fiuctuations of the fiavour 
brane explain the meson spectrum of the field theory. In another interesting 
work this system has been extended to more general settings of Dp-Dg brane 
configurations [13]. 

In the D3-D7 brane set-up the authors of [11] have studied thermaliza- 
tion in the meson sector of the SYM theory. The thermalization is modeled 
by a sudden change in the baryon number. In the gravity dual descrip- 
tion this change is realized by throwing the baryonically-charged fundamen- 
tal strings from the boundary to the bulk. The end points of the strings 
stretched between D3 and D7-branes act as the source for the gauge field on 
the D7-brane. Therefore this set-up provides a time-dependent gauge field 
configuration. Using DBI action this induces a time-dependent metric and 
eventually emergence of an apparent horizon on the D7-brane. Hence the 
baryon injection leads to the horizon formation on the D7-brane which is the 
signal of thermalization in the meson sector of the field theory. 

An intriguing observation of [11] is that the thermalization time in the 
field theory for the D3-D7 system can be written only in terms of a few 
parameters which are: the gauge theory t'Hooft coupling (A = Qym^)) max- 



imum baryon number density (rig) and the inverse of the variation time-scale 
of the baryon number (u). Note that ub and u describe how the baryon 
charge density changes in the system. Therefore the thermahzation time is 
independent of the details of the theory by which we mean the form of the La- 
grangian. The authors generally argue that this possible universal behaviour 
might be generalized to other gauge theories. In this paper we would like to 
examine this idea and generalize their computation to the Dp-Dq system. 

Interestingly for such a general system where the gauge theory is not con- 
formal, the background metric is not AdS and as mentioned above gauge/gravity 
duality can be trusted only in an intermediate range of energy [6] , we observe 
similar universal scaling behaviour. To be more explicit, by universality we 
mean the thermahzation time of the form 

where as we will see in the following sections a and /3 are fixed in terms of p 
using the equation of the apparent horizon. 

The thermahzation time is obtained by studying the dynamics of the 
scalar mesons in Dp-Dg system after the time-dependent baryon injection. 
This can be done by investigating the fluctuations of the transverse and par- 
allel directions of the Dg-brane to the Dp-brane. We will show that although 
the equations of the apparent horizon are different in these two cases, the 
thermahzation time-scales are yet the same. 

We will also observe that even when the background is not supersymmet- 
ric, which means q y^ p, p + 2, p + 4, this behaviour still persists. The only 
change in (1.2) occurs in a and /3. These results confirm, to some extent, the 
claim of [11] that the thermahzation time-scale shows a universal behaviour. 

2 Review on Dp-Dq System 

In the introduction we mentioned that we are interested in studying the 
thermahzation of mesons in strongly coupled field theories in the context of 
gauge/gravity duality. Therefore in this section we will give a brief review 
on Dp-Dq brane configurations that are used to explain the meson spectrum 
in the gauge theory. 

The supergravity solution corresponding to near horizon limit of A^ coin- 



cident Dp-branes is [6] 



ds'^ = H-^{-dt^ + dxD + H^{du^ + u^dnl_p), 
dnl_p = dO"^ + sin^ ednl + cos^ ^rf^r^p^fc, 



(2.1) 



where 



H{u) = {-y~^, e^ = i7^, C,,..., = H-\ (2.2) 

u 

written in the string frame. The Dilaton field is represented by and Coi,..p 
is a (p+ l)-forni field coupled to Dp-branes. The length scale R is defined in 
terms of the string length scale Ig = Vw and the string coupling gg = 6*^°° 

R'-' = W"^ r(^) ggN ll-P. (2.3) 

It is easy to see that for the special case p = 3, (2.1) reduces to AdS^ x 5*^. 
We are considering p < 4 cases because there is no decoupling limit for p > 5 
and hence no dual gauge theory [6]. 

According to gauge/gravity duality, a strongly coupled SU{N) SYM the- 
ory living on the (p+ l)-dimensional worldvolume of N coincident Dp-branes 
is dual to the supergravity on the above background (2.1) in the large A^ limit 
[6]. The isometry group of this background geometry is 50(1, p) x 5*0(9 — p) 
which in the dual gauge theory corresponds to space-time Lorentz symmetry 
and R-symmetry, respectively. The gauge theory coupling constant is related 
to the string coupling and the string length through 

gl^M = 2'Kgg{2'KlsY-\ (2.4) 

One can add a stack of Nf Dg-branes to the above background and study 
it in the probe limit where -^ is small. The low energy effective action for a 
Dg-brane in an arbitrary background is [14] 

S = 5'dbi + 'S'cs , 
5dbi = -Tq I d'i^^i e- V- det(^a6 + B^b + 2'Ka'Fab) , ^2.5) 

5cs = r, /"p[SO(")e^]e2™'^, 



where induced metric Qab and Kalb-Ramond field Bab are given by 

gab = GhlNdaX dbX , 
Bab = BMNdaX dbX 



Mo vN (2.6) 



a,b,... indices are used to describe worldvolume coordinate running over 
0,1,. ..,g. The capital indices M,N,... are used to denote space-time coor- 
dinates. Gmn is the background metric that is introduced in (2.1) and Fab is 
the field strength of the gauge field living on the Dg-brane. Note that in the 
background we consider here the B-field is zero. In the Chern-Simons part, 
(7(") denotes the (n -|- l)-form Ramond-Ramond potential and -P[...] is the 
pull-back of the bulk space-time tensors to the Dg-brane worldvolume. The 

Dg-brane tension is 

1 

In order to embed a probe Dg-brane in the background (2.1) one can 
introduce the following change of coordinates 



-^ = T^^^^TT^m- • (2-7) 



p = M sin 6* , 
a = u cos 9 . 



(2.8) 



Hence the background metric becomes 

ds^ = H-^-dt'^ + dxl) + m (dp^ + p'^dnl + da"^ + a2t/fi?„p_fe) . (2.9) 
This configuration of the Dp-Dg branes can be schematically shown as 

t Xi .. Xd Xd+l .. Xp p Vtk <^ ^7-p-k 

Dp X X X X X XX (2.10) 

Dq X X X X XX 

where 

q = k + d+l. (2.11) 

In order to get a stable background one has to consider the brane configura- 
tions which preserve supersymmetry. This condition dictates that g should 
be p, p+2, p+4 and correspondingly /c is 1, 2 and 3. Note that this embedding 
breaks the initial isometry group of the Dp-brane to 

SO{l,d) X SO{p-d) X S0{k + 1) X S0{8-p-k). (2.12) 



Supersymmetrically embedding Nf flavour Dg-branes into the background 
obtained from tlie near liorizon limit of A^ Dp-branes corresponds to cou- 
pling Nf flavour dynamical matter fields in the fundamental representation 
(dynamical quarks) to the SYM theory which lives on the Dp-brane. This 
reduces the number of supersymmetries of the original theory from 16 to 8. 
The fundamental matter fields which form a A/" = 2 hj^ermultiplet arise as 
the lightest modes of the strings stretched between Dp and Dg-branes. 

The spectrum of the fluctuations of a Dg-brane {Nf = 1) corresponds to 
the mesonic spectrum of a SU{N) SYM theory in p + 1 dimensions which 
is coupled to a hypermultiplet in the fundamental representation [13]. Note 
that in the case where q = p + A, the fundamental fields live on the full p+1 
dimensions of the gauge theory. But for q = p + 2 and q = p they are confined 
to a (c? + l)-dimensional defect where the Dp and Dg-branes coincide. 

3 Thermalization 

The thermalization in field theory is the result of the injection of energy into 
the system. One way to realize it in the context of gauge/gravity duality is to 
inject baryon charge to the strongly coupled field theory. This has been done 
for the first time in [11] where the dynamics of mesons or in other words the 
thermalization of mesons in D3-D7 system has been studied. In this section 
we will elaborate more on this process by generalizing it to Dp-Dg system. 

3.1 Baryon injection 

The injection of baryons into the strongly coupled field theory is presented 
in the gauge/gravity context by throwing byronically charged F-strings on 
the flavour Dg-brane. The end points of the strings stretched between Dg 
and Dp-branes resemble quarks in the gauge theory. The quarks are assumed 
to be massless. They act as a source for the background gauge field on the 
flavour brane. The injection of quarks into the system results in a sudden 
change in the baryon number chemical potential. 

In the AdS/CFT correspondence, static baryon chemical potential is iden- 
tified with the time component of the gauge field in the following way 

du duAt{u) , 
= At{oo)-AtiO) , (3.1) 

8 



where u represents the radial direction in the AdS space^. Therefore in 
order to explain a time-dependent chemical potential which results from the 
injection of quarks to the system we must add the source term 



5', 



current 



j d''^^i^^gAJ\ (3.2) 



to the action (2.5). Note that the currents, J"s, are now time-dependent and 
will produce time-dependent gauge fields. 

Since the end points of the open strings move on the light geodesies, 
without loss of generality, the current components can be fixed to be J^ and 
J*. We assume that the massless quarks are moving along the null line x~ 
defined as 

x± = t±z = tT H'^dp . (3.3) 

Hence the currents are only functions of x~ . 

With the above assumptions the Dg-brane DBI action for the background 
(2.9) in the presence of the source terms reduces to^ 

S = - r,VdVo\{S') f dtdp[H"~'^T"'p''^l - {27ra'yF,l - Af - A^f . 

(3.4) 

To solve the equation of motion for the gauge field we can consistently choose 

f = -H^f = g\x'), (3.5) 

which satisfies the current conservation equation 

VaJ" = 0. (3.6) 

Therefore the classical time-dependent solution for the gauge field is 

i27ia')Ft, = . ^ , (3.7) 

'g^ + i27ra'yp^^H'^'^^^ 



•^Note that this identification is done in the gauge where A^ is zero. 
^Note that Tgj" = ^J"- 



where the following identity has been used 

d^ I p^if-'^l^ , ^'""^^'''*^ I = 0. (3.8) 

It's been already explained that the sudden change in the baryon number 
density can locally model the collision of the heavy nuclei at RHIC The 
baryon number density is defined as n^ = ^ where n is the quark density 
number. In our set-up the quarks live on the {d + l)-dimensional defect 
field theory. Therefore the quark density number is n = y^ where n is 
the quark number (number of open string end points). The time-dependent 
quark number can be calculated by taking the integral of the time component 
of the external source over volume space of the Dg-brane 



n= d'^^y/^r . (3.9) 

Therefore replacing (3.5) in the above equation leads to 

g^x-) = ^^^1^(2^) "^(2™')"^ A nsix-), (3.10) 

where A = Oym-^ ^^^ 7 is a parameter that should be fixed according to the 
convention. For instance in the case of D3-D7 system [11], 7 is ^. 

So far all the calculations have been done for general values of p and q. 
In order to have a stable background in the limit of zero temperature, we 
consider supersymmetric systems where g = p-|-4, p + 2, p and /c = 3, 2, 1, 
respectively. For such configurations we have 

p-g + 2(A;-l) =0 . (3.11) 

Therefore one can see that the p and q dependence in Ftp and g{x~) disap- 
pears and they reduce to 

(2™')i^*p = /,^.l .,,, ' (3-12) 



and 



7 



-(27r)2(2W)^ A 72b(x-), q = p + 4, 



Vol(S3) ' 

9i^') = { voios5)(2^) (2^«')' A nB(x-), q = p + 2, (3.13) 



7 



{2TTa')'^ \ ub^x ), Q = P- 



Vol{S'i) 

10 



Time-dependent gauge field on the Dg-brane modifies the metric observed 
by the open strings. This can be exphcitly seen by using the DBI action. 
We will show in the following section that this field configuration can create 
an apparent horizon on the Dg-brane which signals the thermalization in the 
dual strongly coupled field theory. 

3.2 Apparent Horizon Formation 

Thermalization usually refers to an increase in the temperature of a thermo- 
dynamic system from zero to a nonzero value. We know that in the context 
of gauge/gravity duality the vacuum state in field theory corresponds to pure 
AdS space in the bulk and the thermal mixed state is dual to the AdS-black 
hole. So the horizon formation in the bulk can be considered as a signal 
of thermalization in the dual gauge theory. In order to study the process 
of thermalization in a strongly coupled field theory the gauge/gravity tech- 
niques can be applied [7, 8] . In this framework the horizon formation usually 
refers to the apparent horizon formation since in contrast to the event hori- 
zon it can be obtained locally and the global knowledge of the bulk solution 
is not needed [7]. 

In the Dp-Dg set-up the meson dynamics can be used to study thermal- 
ization in the field theory. The mesons are the open strings (scalar or gauge 
fields) living on the Dg-brane [15]. The Dg-brane can fluctuate in the direc- 
tions which are transverse or parallel to the background Dp-branes. In the 
following subsections we will study the dynamics of these two kinds of scalar 
fluctuations in the presence of a time-dependent gauge field on the Dg-brane. 
It leads to calculating the modified metric on the Dg-brane which is called 
the open string metric [16] . Thus we can obtain the equation of the apparent 
horizon and compute the thermalization time-scale. 

3.2.1 Transverse Fluctuation 

Consider Dg-brane embedding (2.9). The transverse fluctuations are along 
a and flj^p^k directions. Due to the isometry S0{8 —p — k) we assume only 
6a to be nonzero. This plays the role of the scalar meson whose dynamics 
we would like to study. We expand the DBI action (2.5) to quadratic order 
in 6a in the background solution Ftp (3.7). Note that since the RR field does 
not couple to the end point of the string on the brane it will not enter in 
the calculations for the open string metric. Therefore from now on in our 
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calculations we neglect the CS action. Hence the action for the scalar meson 
is 

S = -^ f dtdpd'^x'd'^e'^ v^ g^''da5adb5a, (3. 14) 

where 

. 2-q 



-~gu = rr' [1 - (2W)2Fj]i^i7^t^, 



2-9 P+q-2 



9PP = rr [1 - (2W)^F4]i-^i7^(^, 



p-1 



(3.15) 



~g,, = rr' [1 - {2'Ka'YF']—^H^i^^6,^ 



p+q-2 



-gap = rr [1 - (27ra')'i^J]^p'^^^^^a/3- 



i = l,..,c? is the spatial directions of the id + l)-diniensional defect. ^" 
(a = 1, .., /c) is the angular variable on S^ and ^q,/3 is the metric on the unit 
/i;-sphere. The surface of the presumably apparent horizon, which is defined 
locally as a surface whose area variation vanishes along the null ray which is 
normal to the surface, is given by ^ 



k 

2 
jaoL 



^surface = / d'^x'dH'^ (11^-11 9- 

^ i=\ a=i ' (3.16) 



Tq ..1 / a . /„ ,No o, ^^p-g + 2(fc-l) 



= yrfVol(^*^)^/7t W^2 + i^Tia'Yp^^H' 2 
27ra' '^ 

where this should satisfy 

dVVt=-d. = ^. (3.17) 

Therefore the equation for the apparent horizon reduces to 

^cIe'^P + (p - l)g'^ + z/(2W)2/7"""^§^^p2fc ^ Q^ (3.18) 

where v = ^[{p — 7){p — q + 2k) + ik]. We call this the master equation. 
For supersymmetric cases, the master equation can be simplified further 

4:gg'H'2p + (p - 7)/ + (p - 7 + 2k){2na'Yp^^ = 0. (3.19) 



^One can see that for the supersymmetric cases the apparent horizon surface area 
simphfies to 



T^urfacc = VdYo\iS^)^H-2^g^ + {2T:a'Yp^K 



12 



Since working with z coordinate is more appropriate we replace p with z in 
(3.19) and get 

2(5-p)c/^'2+(p-7)/ + (2W)^(j9-7+2A;)(— ^)^2^i?^^^ =0, (3.20) 

where the length scale R can be written in terms of the t'Hooft coupling 

constant A as 

i?l^ = y2(2W )(27r)^^ r(^^)^ A^. (3.21) 

The formula (3.20) is the main equation that will be used to calculate the 
thermalization time for the transverse fluctuations. 



3.2.2 Parallel Fluctuation 

As it was mentioned earlier in addition to the transverse fluctuations of the 
Dg-brane one can also study the ones which are parallel to the background 
Dp-branes. These parallel fluctuations exist only for g = p + 2 and p. This 
can be easily seen from the conflguration (2.10). Note that for g = p + 2 we 
have d = "p — \ and the only parallel fluctuation is along p direction. For 
q = p there are two directions for parallel fluctuations, p — 1 and p, and the 
fundamental matter lives on the spatial [d = p — 2)-dimensional defect. 
Similar to the transverse fluctuations the DBl action reduces to 

S = -^ f dtdpd'^x'd^e''^g g'^^daSx^'dbSx'' , (3.22) 

where Sx^ represents small fluctuation parallel to Dp-branes. The compo- 
nents of the metric are 

2 — q p — q—4: 

2(9-1) 



-~gu = Tr'[l-i27ra'fFl]—.H 



2 — q p + q — 6 

2(9-1) 



~g,, = Tr'[l-{2na'fF^^]—.H 

~g,, = rr'll - (27rafFj]i^if^(5^5, 



(3.23) 



-gais = rr' [1 - {2'Ka'fFl]^^p'H^^)g^p. 
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Therefore we can now calculate the equation for the apparent horizon area^ 



Kurface = V,Vo\{S^)^H''^\lg^ + (2vr«')V=i/'"'^i^^. (3.24) 

Similar to the previous subsection we apply the condition (3.17) and the 
master equation for the parallel fluctuations reduces to 

2{5-p)gg'z-{p-7)g^-{27va')'^{p-7-2k){^^)^-z^-R''^ =0, (3.25) 

where the length scale R was introduced in (3.21). The master equation 
(3.25) is very similar to the one obtained for transverse fluctuations (3.20) 
except for a couple of sign differences. Note that we have calculated the 
master equation for the supersymmetric cases. 

Before closing this section we would like to emphasize on the observation 
that for the general values of q, not necessarily supersymmetric ones, the only 
change in the equations (3.20) and (3.25) appears in the power of z, in their 
last term. Since p—q+2{k—l) is not zero any more in the non-supersymmetric 
brane configurations the power of z, ^, is replaced by Pz9±J___i^ Note 
that this is a rough argument and we have ignored discussing the stability of 
the backgrounds. We will see from our results in the following section that 
this indicates that the thermalization time scale behaviour does not rely on 
supersymmetry. 



4 Baryon Injection Model 

To be able to determine the time-scale of thermalization, the source function 
g{x^) should be specified. Similar to the function considered in [11], we 
assume 

X" < 

g{x~) = { gmax UJX~ < X- < ^ (4.1) 

9max T^t 



^In supersymmetric cases we have 
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Using the equation (3.10) it shows that the baryon density ub is zero in the 
beginning which resembles the situation before the colhsion of heavy ions. At 
a;~ = it starts changing to reach a maximum value at x~ = - and remains 
constant since then. The maximum value of source function Qmax is given by 
(3.13) where nB{x~) is replaced by its constant maximum value, ub 

9ma. = V^i^(2^)'^(2^"')'"^ ^ ^s- (4-2) 

This is similar to the situation where the baryonically-charged heavy ions 
approach each other and the baryon number chemical potential changes lo- 
cally. 

The time-dependence of the baryon number can be chosen to have a 
general form, {ux~)"', where n > 0^. This is a dimensionless combination 
and guarantees that g{x~) gets its maximal value at a;~ = -. 

We consider the supersymmetric configurations which fall into these cat- 
egories: 

4.1 (p,p + 4) system 

In this set-up of brane configuration we consider q to be p+A where < p < 5. 
As it's been shown in (2.10) there is no parallel fluctuations in this case and 
one needs to consider only the transverse fluctuations. 

4.1.1 Transverse Fluctuation 

If we define g{x^) = Qmax y{x^) where gmax is given by (4.2) and set k to be 
3, the master equation (3.20) becomes 



{p ~ 7)y'^ + 2(5 ~ p)yy' z 



1 , x/5 — n, 12 , , 3(p+i) ^ 7 — », 6 A ^-p 12 l^-"^) 

+ ^(P - l){^)^'i2n)^^ r(^)^— ^^^ = 0. 



^If we assume g{x ) to scale as gmaxi^x )" the change in the thermaUzation time 
appears in the overah power /3 and the power of uj as 



A" 



^th ^ (^ ^) ■ 



IgU 



Note that, in contrast to transverse fluctuations, n — does not give a real solution for 
parallel ones. We will see this in subsection 4.2.2. 
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To warm up, we start with the D4-D8 brane configuration. 
D4-D8 



The master equation is then 

3X^ 



12„,2„2 ^12 



40967r^^-f^n%z 



3y^ + 2zyy' = 0. 



(4.4) 



The location of the apparent horizon on the D8-brane is obtained by solving 
the above equation and determine t in terms of z. We get 



Az 1 



z = 0.16{- 



iriB' 



-;127^12^2^|l^2 


t<z + -, 

UJ 




1 

t> Z+-. 
UJ 



(4.5a) 
(4.5b) 



The curve (4.5a) and the line (4.5b) are plotted in figures 1 and 2 in green 




0.0 0.2 



1.0 1.2 1.4 



Figure 1: The location of the apparent horizon for t < 2; + - in the t — z 
plane. The green curve represents (4.5a) and the blue line shows the tangent 
null ray. The parameters are chosen as a; = 0.5 and ^-^ = -^ for D4-D8 
system. 



and black colours, respectively. The goal of these calculations is to compute 
the thermalization time in the boundary field theory. This in fact is the time 
that a boundary observer starts to see the apparent horizon formation on 
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the D8-brane. For t < z + - the null ray which is tangent to (4.5a) is the 
earliest one that conveys the information to the boundary. It has been shown 
by the blue line in figure 1. The other null rays which cross (4.5a) reach the 
boundary at later times. The tangent point A coordinates are 



ZA = 0.46(^^^ 
So the thermalization time is 



tA = 0.81{- 



X' 



'yuBUJ 



tth = tA + Za^ 



X' 



HbUJ 



(4.6) 



(4.7) 



An interesting observation is that the thermalization time depends only on 
A and the baryon injection parameters ub and u. 




Figure 2: The location of the apparent horizon for t > z -\ — in the t — z 
plane. The black, red and blue fines represent (4.5b), t = z + ^ and tfie 
earliest null ray wfiicfi reacfies tfie boundary. Tfie parameters are cfiosen as 
c; = 3 and (^^)^ = — for D4-D8 system. 



Similar argument can be made for (4.5b). In tfiat case tfie tfiermalization 
time is determined by tfie point A coordinates in figure 2. Point A is wfiere 
tfie line (4.5b) crosses tfie ingoing null ray t = z + -, red line in tfiat figure. 
So tfie null ray wfiicfi passes tfirougfi tfiis point and reacfies tfie boundary 
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gives the earliest time that a boundary observer sees the horizon formation. 
This happens especially if 

(— )^»-, (4.8) 

and thermalization time becomes 

tth = 2zA + - , (4.9) 

where 

1 

tA = ZA + - , (4.10a) 

za = 0.16( )i (4.10b) 

Inequality (4.8) then indicates 

tth~( — )i (4.11) 

Ub 

Dp-D(p+4) 

The above calculation can be generalized to other values of p. We conclude 
that 



<p(^l) <5 



5-p 
2(p-2) \ 2(ll-p) 

«.»~l^l . (4.12) 



2(p-2) \^iJT=p) 

and if ( \^2 J2 I ^ -. we have 



WUUJ^ I CO 



5-p 
2(p-2) \ T^ 

■ A 5-P \ 
t*h ~ ^^ . (4.13) 



n 



B 



p = 1 

Setting p = 1 the system behaves differently from the other values of 
p. This is due to the fact that at p = 1 one of the terms in the master 
equation vanishes and this trivializes the calculations. 



4.2 (p,p + 2) system 

Consider a D(p + 2)-brane supersyminetrically embedded in the background 
obtained from the near horizon geometry of A^ coincident Dp-branes where 
1 < p < 5. k equals 2 and the dual gauge theory is (p + l)-dimensional, 
but the fundamental hypermultiplet has been introduced on a (c? + 1 = 
p)-dimensional defect. There are three and one set of fluctuations along 
transverse and parallel directions, respectively. In the following subsections 
we will discuss them individually. 

4.2.1 Transverse Fluctuations 

For this configuration the master equation (3.20) reduces to 

2{5 -p)yy'z + (p - 7)y'^ 

2{p-3) 8 /A -,A\ 

+ -{p - 3)(27r)^^(^)^ r(^)^ = 0. 

If we repeat the same calculations done in the previous section, the thermal- 
ization time for arbitrary values of p is 



1 <p(^3) <5 



5-p 
2(p-3) \ 2(9-p) 
A 5-P 



5-p 

and if the time scale lies in the limit where ( — ^-A- ^ -, we 



obtain 



i' 



P-5 



( 2(p-3) 
^1 . (4.16) 

p = 3 

Note that similar to p = 1 case in the Dp-D(p + 4) system, choosing 
p = 3 sets one of the terms in the master equation to zero and the 
thermalization time can not be obtained. Interestingly we will see in 
the following that the parallel fluctuation produces the result we expect 
for p = 3. 
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4.2.2 Parallel Fluctuations 

As it was mentioned earlier there is only one parallel fluctuation for Dp- 
D(p + 2) system which is along the x^ direction. The master equation (3.25) 
for k = 2 reduces to 

2{5-p)gg'z-{p-7)g^-{27Ta'f{p-ll){^^)^^z^^R^^ =0 . (4.17) 

Let us start with an explicit example. 
D4-D6 

Setting p = A the apparent horizon equation becomes 



7A^ 
3y' + 2zyy' + ^^^ ^^ = 0. (4.18) 



3g7;-8^2^2 2-8 



This looks similar to the master equation for D4-D8 brane system (4.4). We 
follow the same analogy here. The solution to the above equation which 
specifies the location of the apparent horizon reads 




7A2 1 

t<z + -. (4.19) 



12z^Tr^^'^n%uj 



2 ' 



U 



Note that there is no real solution for z when t > z + -. Thus the thermal- 
ization time is only fixed by the solution (4.19). For t < z+ - the expression 
under the square root must be positive. Fortunately we can always find a 
solution for z that respects this condition. The coordinates of the tangent 
point A to the curve (4.19) will determine the thermalization time. This has 
been fully discussed in D4-D8 case so we state the result. The coordinates 
of the tangent point A are 






(4.20) 



7^/5 UbUJ 

Hence the thermalization time on the boundary reads 



tth = tA + ZAr-i — )^/'. (4.21) 

UbUJ 
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If we set p = 4 in the result for the transverse fluctuation (4.15) we get the 
same thermalization time scale, up to a numeric constant. It is interesting 
that although the apparent horizon equations are different for transverse and 
parallel fluctuations their thermalization times observed on the boundary 
exactly resemble each other. 

Dp-D(p+2) 

The above result can be generalized to other values of p and we get 

One can conclude that the thermalization time for parallel fluctuations and 
transverse ones (4.15) look completely alike for various values of p. An in- 
teresting observation is that the time scale (4.22) works fine for p = 3 case 
in contrast to the transverse one. We can see that for p = 3 there is no 
dependence on the t'Hooft coupling A in the thermalization time scale. 

4.3 (p^p) system 

Consider the brane configuration where q = p- In order to supersymmet- 
rically embed one Dp-brane in the background geometry of A^ Dp-branes 
we choose 2 < p < 5. As it was mentioned earlier the fundamental hyper- 
multiplet is confined to a {d + 1 = p — l)-dimensional surface, k = 1 for 
this set-up. Therefore there are two sets of both, parallel and transverse 
fluctuations which we will elaborate more on them in the following. 

4.3.1 Transverse Fluctuations 

The master equation (3.20) in this case reduces to 

{p - 7)y'^ + 2{5 - p)yy' z 

1 , X.5 — n, 4 ^ , 3(p-3) ,7 — p, 2 A 5-P 4 (4.23) 

+ ^(P - 5)(^)^(27r)^^ r(^)^— ^z^ = 0. 

The thermalization time-scale for general values of p becomes 
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2 <j9< 5 



5-p 

2(P-4) \ 2(7-p) 



ttk - I ^ 1 , (4.24) 



n'^oo^ 



5— p 
2(P-4) \ 2(7-p) 



and in the limit where ( ^ 2 ^2 I ^ -, we have 



W 



p-5 
2(p-4) \ ^ 

^*^ ~ I ^ 1 • (4-25) 



4.3.2 Parallel Fluctuation 

For Dp-Dp system there are two directions for parallel fluctuations. The 
master equation (3.25) reads 

-{p - 7)y^ + 2{5 - p)yy' z 

1 , ,,5-p, 4 ^ 3(p-3) 7-p 2 A^^ 4 (4.26) 

- ^(P - 9)(^)^(27r)^^ r(^)^— ^z^ = 0. 

Interestingly we observe that similar to Dp-D(p + 2) system, the thermal- 
ization time for parallel fluctuations has the same scaling behaviour as the 
transverse fluctuations. Therefore for general p we obtain 



2 <p< 5 



5-p 
2(P-4) \ 2(7-p) 



tth ~ I ^ro^T 1 • (4-27) 



n\uj'^ 



4.4 Concluding Remarks 

We conclude this section mentioning some interesting observations. It's been 
already mentioned that the thermalization time for parallel and transverse 
fluctuations are the same, up to a numeric coefficient, even though the ap- 
parent horizon equations are different. The theories we discussed here are 
not necessarily conformal and their dual background geometries are not AdS. 
But the time-scale behaviour still gets the same form as an AdS background 
with a conformal field theory dual such as D3-D7 system. Even for general 
values of q for which the solution is not supersymmetric the general form of 
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the thermalization time (1.2) persists. This seems to approve the claim that 
the thermahzation time-scale behaves universally. But one may ask to what 
extent this discussion works. We leave it as an open question. 

Let us consider the case where p = 3. As it is well known the dual gauge 
theory is A/" = 4 SYM theory which lives on the 4-dimensional boundary 
of AdS space. If we require to preserve supersymmetry we can add three 
types of probe flavour branes to this background. They are D3, D5, and 
D7-branes. Adding these branes will modify the dual gauge theory to A/" = 4 
coupled to A/" = 2 fundamental hypermultiplet. The dynamical quarks live on 
(1 + 1), (2 + 1) and (3+l)-dimensional defects for g = 3, 5 and 7, respectively. 
Interestingly we see that although the 4-dimensional gauge theory is the same 
for all of these brane set-ups these theories produce different time scales for 
thermalization which are 

(;^)^ D3-D7, 

{^)K D3-D5, (4.28) 

)i D3-D3. 



tth 



r^gu'^X' 



Therefore we can conclude that the thermalization time scales differently 
with the t'Hooft coupling, A. The main difference between these theories 
come from the fact that the dynamical quarks are confined to defects with 
different dimensions. 

Another intriguing observation is when we consider the dimension of the 
defect to be (2 + 1). The brane configurations to get it are D2-D6, D3-D5 
and D4-D4. One can see that in all of these cases the dependence on A in the 
thermalization time disappears. It is very interesting that for this specific 
choice all the brane configurations give the same value for the thermalization 
time which is ( 2^ 2 )^. Note that this doesn't happen in the other dimensions 
for the defect such as (1 + 1) and (0 + 1) dimensions. 

In this paper we have studied thermalization for the scalar mesons. One 
can generalize this calculation and study the dynamics of the vector mesons. 
Moreover the universal behaviour of the thermalization time-scale can be 
investigated in other gravity backgrounds such as anisotropic Lifshitz ones. 
It is also interesting to study the time-scale of the change in the temperature 
if one starts from thermal field theory, dual to black hole background, instead 
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of zero-temperature one and then injects energy into it. We postpone these 
to the future works. 
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